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The  direct-interaction  approximation  for  turbulence  is  extended  to 
predict  the  covariance  and  average  Green’s  function  of  a  generalized  velocity 
u(x,t|r).  The  latter  is  defined  as  the  velocity  measured  at  time  r  in  the 
fluid  element  which  passes  through  the  ooint  x  at  tire  t.  The  resulting 
formulas  for  triple  moments  involve  integrals  over  the  Eulerian  time-history 
of  the  fluid.  The  approximation  is  then  altered  so  that  the  integrals  are 
instead  over  Lagrangian  histories,  measured  aionr  the  panicle  paths.  The 
alteration  is  made  necessary  and  is  uniquely  determined  by  reauiring  simul¬ 
taneously  the  consistency  properties  that  energy  be  conserved;  that  there 
exist  formal  inviscid  equioarrition  solutions;  and  that  the  dynamics  exhibit 
invariance  under  a  class  of  random  Galilean  transformations.  In  the  altered 
approximation,  the  relaxation  times  associated  with  energy- transfer  are 
Lagrangian  memory  times  determined  by  the  viscous  and  pressure  forces.  As 
a  result,  the  approximation  yields  the  Kolmogorov  inertial-  and  dissipation- 
range  laws.  The  corresponding  approximation  for  convection  of  a  passive 
scalar  field  yields  seme  exact  results  of  Tavlcr  «:nd  yields  Richardson's  law 
for  the  relative  diffusion  of  two  particles. 


1.  INTRODUCTION 


It  was  stressed  in  an  earlier  paper^  that  Euierian  moments  do  not  provide 
an  appropriate  description  of  the  convection  of  small  spatial  scales  of  turb- 
ulence  by  large  spatial  scales.  Knowledge  of  the  Euierian  velocity  covariance 
alone  does  not  permit  discrimination  between  two  importantly  different 
situations.  In  one,  the  time-dependence  of  the  small-scale  features,  measured 
at  fixed  points  in  space,  is  due  to  internal  distortion.  In  the  other,  the 
time-dependence  is  due  to  the  small  scales  being  s*ept  along,  almost  undis¬ 
torted,  by  the  large-scale  motion.  In  the  first  situation,  the  measured 
characteristic  times  of  variation  are  relevant  to  the  transfer  of  energy 
among  the  small  scales,  while  in  the  second  this  is  not  so.  This  suggests 
that  closure  approximations  which  involve  only  low-order  Euierian  moments  do 
not  retain  sufficient  information  to  represent  properly  the  energy  transfer 
among  small  scales  which  are  connected  by  large  scales.  In  Ref.  1  it  was 

shewn  in  detail  that  this  difficulty  makes  it  impossible  for  the  direct- 
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interaction  approximation  *  and  a  related,  higher  Euierian  approximation 
to  predict  correct  inertial  range  dynamics. 

A  generalized  velocity  u(x,tjr)  may  be  defined  as  the  velocity  measured 
at  time  r  within  the  fluid  element  which  passes  through  the  point  x  at  time  t. 
The  function  u(x,t|t)  coincides  with  the  Euierian  field  u(x,t),  while  u(x,t|r) 
considered  as  a  function  of  r  at  fixed  x  and  t  is  equivalent  to  the  usually 
defined  Lagrangian  velocity.  The  two  situations  posed  in  the  preceding  para¬ 
graph  can  give  similar  t  dependence  of  u(x,t)  at  fixed  x.  However,  internal 
distortion  of  the  small  scales  implies  variation  of  u(x,t|r)  with  both  t  and 
r,  while  mdistorted  convection  of  these  scales  implies  that  the  convection 
induces  a  change  of  the  t  dependence  but  no  change  in  the  r  dependence  at 
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fixed  x.  This  suggests  that  knowledge  of  the  covariance  of  the  full  function 
u(x,t$r)  say  permit  sufficient  discrimination  between  the  two  cases. 

In  the  present  paper,  an  alteration  of  the  direct-interaction  approximation 
is  formulated  which  £ives  closed  statistical  equations  involving  the  covariance 
of  u(xvtir)  and  the  average  response  of  u(x,t|r)  to  infinitesimal  perturbations. 
Construction  of  the  final  approximation  involves*  two  stages.  First,  the  direct- 
interaction  approximation  for  u(x,t  J  r)  is  constructed  by  a  straightforward 
extension  of  previously  used  techniques.  This  yields  expressions  for  the 
triple  moments  in  terms  of  the  covariance  function  and  the  average  response 
function.  At  this  stage,  the  equation  for  evolution  of  the  Eulerian  covariance 
is  exactly  the  same  as  in  the  purely  Eulerian  direct -interact ion  scheme. 

The  equations  so  formed  preserve  certain  fundamental  properties  of  the 
exact  dynamics:  Conservation  of  energy  by  the  nonlinear  interaction,  maintain- 
ance  of  the  incompressibility  of  the  Eulerian  field,  invariance  of  total 
stress-energy  under  the  transformation  from  Eulerian  to  Lagrangian  coordinates, 
and  the  existence  of  formal  inviscid  equipartition  equilibrium  states.  The 
inability  of  the  direct-interact loo  equations  correctly  to  represent  convection 
effects  is  displayed  in  sharp  form  by  failure  to  preserve  a  further  fundamental 
property  of  toe  exact  dynamics:  Invariance  under  random  Galilean  transformation. 
Suppose  that  each  flow  in  the  statistical  ensemble  is  subjected  to  a  transla¬ 
tional  velocity  which  is  constant  in  space  and  tine  but  which  has  a  Gaussian 
enoeabl*  distribution.  Clearly  the  internal  dynamics  of  the  turbulence  is 
unaffected,  and  this  is  expressed  analytically  by  certain  transformation  laws 
for  the  statistical  functions  under  the  random  translation.  These  laws  are 
badly  violated  by  the  direct-interaction  equations. 
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The  second  stage  in  the  procedure  is  to  alter  the  direct-interaction 
expressions  for  the  triple  moments  in  such  a  way  ?s  to  keep  the  conservation, 
incompressibility ,  invariance,  and  equilibrium  properties  already  incorporated 
and  to  realize,  simultaneously,  a  restricted  form  of  invariance  to  random 
Galilean  transformation.  The  unaltered  approximation  for  the  triple  moments 
involves  space-time  integrals  which  excess  memory  and  relaxation  effects  in 
the  turbulent  motion.  The  integrals  have  the  form  of  Euierian  time-histories 
(histories  at  fixed  points  in  space  which  then  are  integrated  over  space).  In 
the  eltered  approximation,  the  space-time  integra’-i  are  changed  so  that  they 
instead  are  over  Lagrangian  histories  (fluid-element  space-time  trajectories). 
The  effective  relaxation  or  memory  times  are  now  measured  in  coordinates 
moving  with  the  flow  instead  of  in  fixed  coordinates,  and  the  evolution  of 
the  Euierian  velocity  covariance  is  now  inextricably  coupled  to  that  of  the 
full  Lagrangian  covariance. 

The  prescription  for  the  alteration  is  heuristic.  However,  it  involves 
no  arbitrary  constants  or  functions,  and  seems  uniquely  determined  by  the 
required  invariance,  conservation,  and  equilibrium  properties.  One  important 
property  of  the  original  approximation  is  not  shown  tc  survive  the  alteration. 
The  direct- interact ion  equations  are  exact  statistical  equations  for  a  certain 
dynamical  model.  This  guarantees  that  the  results  for  the  covariances  are 
realizable  as  averages  over  some  ensemble  of  velocity  fields  and  that,  conse¬ 
quently,  gross  unphysical  predictions  such  as  negative  spectra  are  precluded. 
The  guarantee  does  not  presently  exist  for  the  altered  approximation.  The 
solutions  so  far  found  give  no  hint  of  unphysical  behavior  and,  in  every  case, 
represent  improvements  over  the  predictions  of  the  unaltered  equations. 
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The  new  equations  are  called  the  LHDI  (Lagrangian-History  Direct- 
Interaction)  approximation.  They  are  constructed  first  for  the  convection  of 
a  passive  scalar  field  by  a  prescribed  turbulent  flow,  and  then  for  the 
dynamics  of  the  turbulent  velocity  itself. 

The  UiDI  equations  appear  to  yield  inertial  and  dissipation  ranges  at 
high  Reynolds  numbers  which  obey  Kolmogorov's  laws.  This  is  because  the 
effective  relaxation  tines  in  the  energy- transfer  function  are  Lagrangian 
memory  times.  The  approximation  taay  in  fact  be  regarded  as  a  mathematically 
precise  (but  dynamically  approximate)  embodiment  of  Kolmogorov's  original 
idea  that  the  dynamics  of  straining  and  energy-transfer  should  be  examined  in 
coordinates  which  move  locally  with  the  fluid.  In  the  case  of  the  scalar 
field,  the  LHDI  equations  yield  some  exact  results  of  Taylor  for  the  dispersion 
of  particles  by  homogeneous,  stationary  turbulence.  When  the  velocity  field 
has  a  Kolmogorov  inertial  range,  the  scalar  equations  yield  Richardson's  law 
for  the  relative  diffusion  of  txo  particles. 


2.  GENERALIZED  VELOCITY  FIELD 

Let  u(x,t)  be  the  velocity  at  time  t  at  the  point  x  in  a  fixed  Cartesian 
coordinate  system  (Eulerian  velocity).  Define  the  field  u(x,t|r)  by 

u(x,r]r)  *  u(x,r),  (2.1) 

(>/?t  ♦  u( x ,t ) * V}u( x , t | r )  =  0.  (2.2) 

Eq.  (2.2)  bolds  for  all  t  and  r  (t  >  r  and  t  <  r).  No  spatial  boundary  con¬ 
ditions  on  u(x,t:r)  are  needed,  or  may  be  imposed,  for  t  i  r.  Eqs.  (2.1)  and 
(2.2)  imply  that 


(2.3) 


u(x,t | r )  =  u(x-(,(x,t  |r),r), 

where  £(x,tjr)  is  defined  as  the  displacement,  during  the  interval  r  to  t,  of 
the  particle  which  arrives  at  (x,t).  Thus,  u(x,t|r)  is  the  velocity  at  tine 
r  of  the  particle  which  arrives  at  x  at  tine  t.  In  what  follows,  u(x,tjr) 
will  be  called  the  generalized  velocity.  The  time  argument  preceding  the 
vertical  bar  will  be  called  the  labeling  tine  and  that  following  the  bar  will 
be  called  the  measuring  tine.  If  (x,t)  is  taken  as  some  particular  point  (a,tQ) 
and  t.  is  an  initial  tine,  then 

u(a,tQjr)  =  w(a,r),  (2.4) 

4 

where  w(a,t)  is  the  Lagrangian  velocity  as  usually  defined.  Considered  as  a 
function  of  t  at  fixed  r,  u(x,t|r)  gives  the  velocities  neasured  at  time  r  of 
the  fluid  elements  which  pass  through  a  given  point  £  at  various  times  t. 

How  let  u(x,t)  obey  the  incompressible  Havier-Stokes  equation 

[3/at  -  vV2  u(x,t)  *Jju(x,t)  =  -?p,  (2.5) 

V*u(x,t)  =0,  (2.6) 

where v  is  kinematic  viscosity  and  p  is  kinematic  pressure.  Then  (2.1),  (2.2), 
(2.5),  and  (2.6)  form  a  complete  set  which  determine  u(x,t)  end  u(x,t!r)  when 
the  spatial  boundary  conditions  and  initial  conditions  on  u(x,t)  are  specified. 
Fig.  1  illustrates  the  path  which  u(x,t!r)  follows  in  the  (t,r)  plane  as  it 
evolves  from  a  specified  initial  field  u(x,tQit0).  The  evolution  along  the 
diagonal  from  (t0,tQ)  to  (r,r)  obeys  (2.5)  and  (2.6),  and  the  evolution  from 
(r,r)  to  (t,r)  obeys  (2.2). 

The  invariance  oroperties 

-r—  Ju.  (x,t  ir)u.(x,t  lr)d  x  =  0,  etc.  (2.7) 

3t  '  l  —  1  ~  * 


|f  /u.(x,t|r)d3x  =  0, 
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C.  1  i  — .  — —  j;  i  ..  £ _ ^  ^  a 

iviiw  ttauny  iiw  \  t  *  t  )  auu 


(2.6),  provided  that  the  normal  cotnoonent  of 


the  Eulerian  velocity  vanishes  on  the  boundaries.  Eq.  (2.7)  reflects  the 
fact  that  u(x,t|r)  is  a  relabeling  of  the  values  of  the  Eulerian  field  u(x,r) 
according  to  a  coordinate  transformation  defined  by  the  displacement  E(x,tir). 
In  the  case  of  an  incompressible  Eulerian  field,  the  Jacobian  of  the  trans- 

4 

formation  equals  one.  In  general. 


v»u(x,t|r)  i  0  ('•  i  r).  (2.8) 

«•**■  *** 

It  is  easily  seen  from  (2.2)  that  (2.6)  does  not  imply  a  divergence less 
generalized  velocity. 

Consider  next  the  passive  scalar  field  *(x,t)  which  satisfies 

(3/it  -  «V2)y(*,t)  =  -  u(x,t jt)«7y(x,t) ,  (2.9) 

where  xis  the  kinematic  diffusivity  and  (2.6)  is  assumed.  A  generalized  field 
♦(x,t|r)  may  be  defined  by 

y(x,r|r)  =  P(x,r),  (2.10) 

3*(x,tjr)/3t  *  -  u(x,t |t)-V*(x,t|r).  (2.11) 

It  satisfies  invariance  relations  similar  to  (2.7).  If  t  :  o,  it  is  clear 
from  (2.9)-(2.11)  that  v(x,tjr)  is  independent  of  r,  which  expresses  the 
constancy  of  the  scalar  concentration  along  the  particle  trajectories. 


3.  RESPONSE  TO  PERTURBATIONS 

Construction  of  the  direct-interaction  approximation  requires  the  intro¬ 
duction  of  the  Green's  functions,  or  response  functions,  which  describe  the 
propagation  of  arbitrary  infinitesimal  perturbations  of  the  system.  C are  is 
needed  in  order  to  define  these  functions  in  a  consistent  fashion  for  the 
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generalized  fields.  Consider  first  the  scalar  field.  The  Lulerian  Green's 


function  may  be  defined  by 


G(x,t  ;x'  ,t' )  =  6*(x,t)/6f(x' ,t' ), 
G(x,t;x',t')  =0  (t  <  t'). 


(3.1) 


where  f(x#t)  is  an  arbitrary  source  term  added  to  the  rignt-hand-side  of 
(2.9)  and  6/6  denotes  functional  differentiation.  It  follows  that 


(3/3t  -  <^)G(x,t;x*  ,t’)  =  -u(x,t|t)*VxG(x,t;x’ ,t')  (t  >  t'),  (3.2) 


G(x,t  * ;x' ,t' )  =  6  (x-x' ) , 


(3.3) 


where  6  (x-x' )  is  the  three-dimensional  Dirac  function. 

The  Green's  function  for  the  generalized  field  may  be  defined  by 

G(x,t|r;x' ,t' |r' )  =  5v(x,t | r)/5f (x' ,t' 1 r' ) 
G(x,t|r;x' ,t' jr')  =0  (r  <  r'). 


(3.4) 


where  f(x,t|r)  is  an  arbitrary  source  term  added  to  the  right-hand  side  of 
(2.11).  Eq.  (3.4)  is  to  be  interpreted  in  the  following  way:  The  perturbation 
in  t  propagates  from  (t',r')  to  (r',r' )  according  to  the  equations 


3G(x,t jr' ;x' ,t* jr' )/3t  =  -  u(x,t it)»VxG(x,t jr' ;x* ,t ' jr' ) , 
G(x,t* jr' ;x' ,t* |r* )  =  63(x-x' ). 


(3.5) 

(3.6) 


Then  [note  (2.10)]  it  propagates  from  (r' #r')  to  (r,r)  according  to 

(3/3t  -  «v^)G(xtt jt;x’ ,t* |r* )  =  -u(x,t jt)»VxG(x,t |t;x* ,t' jr* ).  (3.7) 

Finally,  it  propagate",  from  (r,r)  to  (t,r)  according  to 


JG(x,t jr;x' ,t' |r' )/3t  =  -u(x,t)t)«7  G(x,t |r;x' ,t' |r* ). 


(3.8) 


Eq.  (3.1)  say  be  considered  a  special  case  of  (3.4); 


G^.tjx*  ,t')  =  G(x,t  |t;x' ,t' }t*  ). 


(3.9) 
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Also, 

G(x,tjt;x* ,t' |r' )  =  6*(x,t )/6f(x' ,t ' jr' ) , 

G(x,t|r;x'  ,t*  jt' )  =  6*(j^,t|r)/6f(xf,tJ).  (3.10) 

Consistency  requires  that  f(x,t|r)  be  taken  as  a  source  tern  in  (2.11)  when 
t  4  r  and  in  (2.9)  when  t  =  r. 

The  generalized  Green's  function  has  the  following  physical  significance. 
Suppose  that  a  perturbation  is  externally  imposed  at  tine  r*  in  the  scalar 
concentration  at  the  fluid  elenent  which  passes  through  point  x*  at  tine  t'. 

A  .  . 

Then  G(x,t jr;x' ,t' |r* )  gives  the  resulting  perturbation,  at  tine-of-neasurenent 
r,  in  the  scalar  cooentration  at  the  fluid  elenent  which  passes  through  point 
x  at  tine  t.  Fig.  2  illustrates  the  way  in  which  the  perturbation  propagates 
in  the  (t,r)  plane.  Note  that  there  are  no  restrictions  on  the  labeling  tines 
t  and  t '  but  that  propagation  is  always  in  the  direction  of  increasing  tine 
of  neaaurenent  (r  >  r'). 


The  velocity- field  Green's  function  tensor  nay  be  defined  by 
G,  .(x,t  |r;x* ,t* |r*)  =  6u.(x,t jr)/6f  (x' ,t' |r' ) , 
G..(x,tjr;x*,t' jr')  =0  (r  <  r’). 


(3.11) 


Here  f(x,t|r)  is  a  force  tern  added  to  the  right-hand  side  of  the  equation 
of  notion  for  u(x,tjr),  if  t  ^  r,  and  to  the  right-hand  side  of  the  equation 
of  notion  for  u(x,t)  *  u(x,t|t)  if  t  =  r.  The  path  of  propagation  of  the 
perturbation  is  the  sane  as  for  the  scalar  case.  In  contrast  to  the  scalar 


case,  the  equation  of  notion  for  the  velocity  field  is  nonlinear.  Therefore, 
f(x,t)r)  nust  be  taken  as  an  infinitesinal  if  the  result  for  G^..  is  to  be 
independent  of  f . 


Special  consideration  nust  be  given  to  inconpressil ility.  In  an  in- 
conpressible  fluid  not  all  initial  perturbations  of  the  velocity  field  are 
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poesible.  If  only  the  Eulerian  field  is  considered,  the  restriction  that 
5f(x*  ,t')  be  divergence  less  can  be  i  noosed.  However,  u(x,t|r)  is  not  diver¬ 
genceless  in  genaral.  The  possible  disturbances  6u(x* ,t'|r')  are  those  which 
propagate  to  the  diagonal  of  tne  (t,r)  olane  to  give  a  divergence less  Eulerian 
field,  and  this  does  not  correspond  to  a  sisrole  restriction  on  6f(x' .t* lr* ). 

This  situation  nay  be  handled  correctly  by  Beans  of  a  formal  artifice.  Let  a 
fictitious  curlfree  part  of  the  Eulerian  field  be  admitted  and  write 

C  p 

u(x,t|r)  =  u  (x,tjr)  ♦  u  (x,t!r)  (all  t  and  r),  (3.12) 

where 

V-u^U.tj  r)  =  0,  V  *  uC(x,tjr)  =  0.  (3.13) 

Eqs.  (3.12)  and  (3.13)  imply 

u^(x.tjr)  =  ?. .(y)u.(xyt  jr),  u^(x,t|r)  =  lK^(v)u .(x,t |r),  (3.14) 

where 

Pi.(y)  =  6..  -  n..(v),  (3.15) 

and,  for  any  g(x) , 

!!..(V)g(x)  =  32/3x  3x.  |D(x,y)g(y)d3y.  (3.16) 

ij—  i  i  4  —  —  — 

The  integration  in  (3.16)  extends  over  the  whole  volume  occupied  by  the  fluid. 


D(x,y)  has  zero  normal  derivative  on  the  boundaries  and  satisfies 

y^D(x,y)  =  63(x-y).  (3.17) 

How  replace  (2.2)  and  (2.5)  by 

i»u(x,t  |r)/3t  =  -uS(x,t  |t)*Vu(x,t  |r)  ,  (3.18) 

(3/3t  -  v?~)uS(x,t |t)  =  -uS(x,t jt)»7uS(x,t |t)  -  7p,  (3.19) 

auC(x,t|t)/at  =  0.  (3.20) 


Clearly  in  these  equations  the  fictitious  field  u^ix.tjt)  does  not  affect  the 

s  sc 

evolution  of  u  (x,t|t)  or  the  parts  of  u  (x,tjr)  and  u^Xjtjr)  induced  by 
S  C 

u  (x,t|t).  If  u  (x,t|t)  is  initially  zero,  it  stays  zero.  Therefore  no 


i 


violence  has  been  done  to  the  physics.  However,  (3.20)  now  defines  the 
propagation  of  ''impossible"  disturbances  on  the  diagonal  of  the  (t,r)  plane, 

A 

so  that  the  function  i r;x' ,t * | r' )  is  well  defined  for  arbitrary  initial 

perturbations.  The  "impossible"  perturbations  actually  are  of  physical 
interest.  Their  propagation  provides  a  measure  of  the  distortion,  relative 
to  a  Cartesian  system,  of  a  coordinate  system  which  moves  with  the  fluid.  In 
particular,  the  perturbation  6u  (x,t*|t*)  induced  by  a  divergence less  force 
«A5.fi  r')  provides  such  a  measure. 

It  is  convenient  now  to  eliminate  the  oressurv  from  the  equations.  When 
there  is  no  shear  at  the  boundaries  (e.g.,  boundaries  at  infinity)  the  elimi¬ 
nation  gives 

0/3t  -  W2)u.(x,t|t)  =  -  —  P. .  (V)[u^(x,tit)u5(x,tit)J,  (3.21) 

1  —  <C  lj®  ~  j  ~  B  — ' 


where 

P..  (V)  =  P..(V)3/3x  ♦  P.  (V)3/3x. .  (3.22) 

13®  —  lj  —  m  xb  —  3 

Hote  that  (3.21)  is  equivalent  to  (3.19)  and  (3.20)  together.  This  equation 
will  be  adopted  hereafter.  More  general  boundary  conditions  can  be  handled 
by  the  methods  of  Ref.  3. 

The  eauations  for  G..  can  now  be  obtained  bv  introducing  force  terms  on 

13 

the  right-hand  sides  of  (3.18)  and  (3.21)  and  performing  the  functional 


differentiation.  The-  results  are 
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where 


G^  .(x,t  |  r;x'  ,t '  |  r*  )  =  P.  (V  )G  .(x,t  | r;x'  ,t*  I  r* ) , 
ij  -  ~  ltn  -x  mj  -  ~ 


(3c  26) 


The  property  P. .  (?)  =  P.  .(?)  is  used  in  writing  (3.25). 
K  J  ljm  "  lnsj  ~ 


4.  STATISTICAL  FUNCTIONS 

Consider  a  distribution  of  the  fields  over  statistical  ensemble  and 
assume  that  the  mean  fields  vanish  for  all  argument  values.  The  simplest 
statistical  quantities  are  then  the  covariances 

*(x,t  jr;x'  ,t ' jr' )  =  (»(x_,t!r)Kx'  .t'fr’)), 

U. .(x,t |r;x'  ,t* |r* )  =  <u  (x,t  |r)u  (x*  ,f  |r’  )>  ,  (4.1) 

X  J  A  J 


which  have  the  symmetry  properties 

»(x,t |r;x* ,t* |r' )  =  T(x' ,t* |r* ;x,t |r), 
U.-Cx.tlrjx' ,t' jr’ )  =  U_..(x',t'|r’;x,t|r). 

Ensemble-averaged  Green's  functions  may  be  defined  by 
G(x,t|r;x\t'  |r')  =  <G(x,t  |r;x»  ,t’ jr' )>, 
G.^(x,t  |r;x*  ,t'  |r' )  =  (g.  .(x,t |r;x'  ,t’  |r* )). 


(4.2) 


(4.3) 


They  satisfy 


G(x,t* |r* ;x' ,t' lr' )  =  4  (x-x'}, 

G. ,(x,t' |rf ;x' ,t' ir’ )  =  6,.53(x-x'). 

~  ii  — - 


(4.4) 


In  correspondence  to  (3.14),  U..  may  be  decocoosed  into  the  tensors  U.. 

’  11  7  IT 


MC  ,.5  C  nSS  „SC  nCS  „CC  ,  ,  K 

U...  U...  U..t  U...  U .  . ,  'J...  U..  defxnec  by 
ii'  xi*  i  j  *  ir  it  *  n*  n  7 


U. .(x,t ir;x’ ,t' lr' )  =  P.  (?  )U  .(x,t !r;x' ,t' Ir' ) , 
U.^(x#t |r;x* ,t* lr' )  =  H.  (?  )U.  (x,t lr;x' ,t * ir' ) , 

ij  ~  -  *  1B1  X  Iff  - 
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The  notation  used  is  that  a  superscript  S  or  C  refers  to  the  index  directly 
beneath  it.  The  corresponding  decomposition  of  G.^  is 


6ij<5tt!r;*\f|r')  =  P«.C2w><v4U.t  {r;x*  ,f  |r« >, 


xss  -x  si 


t‘!r»)  =  G  U,tjr;x’,t'jr')n  (7  ),  etc.  (4.6) 

*  i  a*°  n  j  x 

The  follow isg  example  illustrates  the  use  of  (4.6).  The  vector 


/sljtx.tjrjx* ,t* jr* )5u.(x* ,t’ ir* )d3x* 

gives  the  mean  response(5u  (x,tir))>  of  the  shear  field  when  the  imposed 
(statistically  sharp)  disturbance  produces  an  initial  Derturbation  fu(x’ ,t * jr' ). 
The  vector 

/Gi^(x,t ir;x* ,t* (r* )iu.(x' ,t' |r* )d3x* 

gives  the  total  response  to  the  covoressive  part  5uC(x*  ,t’  ir’ ) 

of  the  initial  perturbation. 

When  the  ensemble  is  homogeneous,  the  covariances  and  averaced  Green's 
functions  depend  on  the  arguments  x  and  x*  only  in  the  cosabination  x-x! .  In 
this  case  the  wave vector  functions 


fUjtlritMr')  =  (2v)”3/T(x,t  ir;x'  ,t '  |r'  )e-1-  )d(x-x' ) , 

G(k;tir;t'ir')  =  /c(x,t  jr;x*  *t '  |r*  )e'^*(-"-  ^(x-x' ) , 

U_(k;t|r;t»  |r')  =  (2»>"3Ju. .(x#t  |r;x'  ,t '  |r*  )e”1*,<*'~' ^(x-x' ) , 
GijUjtlr-.t’lr')  =  /G_(x,t|r;x'  ,t'  ir' )e'1-*(~"-  ^(x-x' )  (4.7) 

form  a  natural  descriotion.  The  normalization  chosen  gives 
/T(^;t|r;t' |r')d3k  =  f(x,t !r;x,t* Ir* ) , 

/u.^(k;t|r;t' |r’ )d3k  =  U..(x,t !r;x,t' jr* ) 


(4.8) 
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and 

G(k  ;t  *  j  r*  ;t  *  |  r’  )  =  1,  G .  { k ;t 1 \ r' ;t ’ | r* )  =  5...  (4.9) 

When  the  distribution  has  refiectional  symmetry,  the  functions  are 
invariant  to  interchange  of  %  and  .  Then  it  follows  that 


f,SC  ,,CS  „SC  CS 

U .  .  =  U .  .  =  0 ,  G .  .  =  G .  .  =  0 . 
13  13  13  i: 


(4.10) 


for  all  argument  vaiues.  In  a  homogeneous  distribution  with  refiectional 
symmetry,  the  decomposition  of  U..  is  completely  described  by 


n 


U^.(k;t jr;t' |r' )  =  P.  (k)U  (k;t |r;t* |r» ) , 
13  -  ic  ®3 

u9.(k;t  |r;t'  |r' )  =  H.  (k)U  .(k;t 1 r;t ' jr* ) , 

13  ~  ic  ir>j  ~  ’  * 


where 


P..(k)  =  6..  -  II.. (k),  H..(k)  =  k.k.k 

13  -  13  13  13  ~  13 


-2 


(4.11) 

(•*.12) 


The  decomposition  of  G. .  is  similar. 


5.  DIRECT-INTERACTION  EQUATIONS :  SCALAR  FIELD 

Assume  tnat  the  fields  are  normally  distributed  at  an  initial  time  tQ, 
with  zero  means  and  zero  correlation  between  scalar  and  velocity  fields.  Th«*n 
the  initial  distribution  is  conoietelv  specified  by  t(x,t.Jt_ ;x' ,t_  It  )  and 
U„(x,t0  lt0;x’  ,tQ  |tQ).  Assume  that  the  r>~an  fields  stay  zero  for  ail  tines. 
The  relations  in  Sec.  3  lead  to  the  equations  of  motion 


(3/3t 

-  *^)t(xft  | t;x '  ,t'  |r' ) 

=  S(x,t 

1 1 ;x* ,t ' 

ir’). 

(5.1) 

O/at 

-  «7^)G(x,t lt;x’ ,t’ jr* ) 

=  rl(x,t 

It ;x* ,t  * 

1  r' )  , 

(5.2) 

3?(x,t |r;x' ,t*  |r' )/3t 

=  S(x,t 

i  r;x' ,t  * 

ir’). 

(5.3) 

3G(x,t|r;xs,t’ |r’)/3t 

=  n(x,t 

lr;x\f 

Ir’), 

(5.4) 
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where 

S(x,t jr;x* ,t* jr* )  *  -^u(x,t |t)«Vxp(x,t jr)*(x* ,t* jr' )),  (5.5) 

H(x,t jrjx* ,t* |r* )  =  -{u(x,t |t)«?xG(x,t |r;x* ,t* jr’  )) .  (5.6) 

Tbe  direct -interaction  procedure  yields  approximations  for  S  and  H  in 

tents  of  f,  G,  and  U..  and  thereby  produces  a  closed  set  of  equations.  The 

2  3  5  6 

method  has  been  described  for  £ulerian  fields  in  several  papers.  '  *  ' 
extension  to  the  generalized  fields  can  be  made  straightforwardly  by  using 
tbe  formalism  developed  in  the  preceding  Sections.  The  simplest  way  to  con¬ 
struct  the  approximation  is  an  algorithm  based  on  iteration  expansion  of  S 
and  H.  Consider  the  zeroth-order  equations  obtained  by  setting  equal  to  zero 
the  right-hand  sides  of  (2.9),  (2.11),  (3.2),  (3.5),  (3.7),  and  (3.3).  The 
solutions  of  these  equations  are 

G(0)(x,t|r;x',t*|r»)  *  [4«(r-r*  )r3/2exp[-|x-x»  |2/4c(r-r*  )], 

*t0)(x,t|r)  *  /G(0)(x,t|r;x',t0!t0Mx',t0it0)d3x*  (5.7) 

Ain)  MO)  (0)  (0)  i 

C»  displays  no  statistical  fluctuation  (G  =  G  ).  Since  *  (x,t|r) 

is  a  linear  functional  of  the  initial  field,  it  is  normally  distributed.  Note 
that  G^°\x,t  |r;x*  ,t*  |r' )  and  f^°^(x,tlr)  have  values  independent  of  the 
labeling  times  t  and  t'.  The  similar  procedure  for  the  velocity  field,  setting 
the  right-hand  sides  of  (3.21),  (3.24),  and  (3.25)  equal  to  zero,  gives 

£^>(*,t|r;x,,t,|r,>  r  6_0»v(r-r'  )j"3/2expC-|x-x'  !2/4v(r-r» )], 

u!°^(x,t!r)  *  /G^?^(x,t |r;x* ,t  ltn)u.(x*,t  It  )d3x*.  (5.8) 

X  ^  1]  ^  0  0  j  uu 

/\  A 

formally  exact  expansions  of  G,  *,  G.. ,  and  u  as  functional  power  series 
in  the  zeroth-order  functions  can  now  be  generated  by  introducing  the  actual 
right-hand  sides  of  the  equations  of  motion  as  perturbations  and  iterating. 
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Tnese  series  may  be  substituted  into  (5.5)  and  (5.6)  and  the  averages  evalu¬ 
ated  Dy  tne  rule  for  moments  of  a  normal  distributi on.  The  results  are  in¬ 
finite  functional  power  series  which  express  S  and  ri  in  tents  of  g!?^, 

and  Clearly  only  pur-.  Lulerian  velocity  functions  appear  in  the 

results  since  only  the  Eulerian  velocity  enters  the  equations  of  motion  for 
v(x,t |r) 

Tne  lowest--,rder  terms  in  the  series  for  S  and  H  are  respectively  linear 
and  bilinear  in  G^°^;  they  do  not  involve  G^?\  The  direct-interaction 
approximation  is  constructed  by  retaining  only  the  lowest-order  terms  and,  in 
them,  replacing  all  zeroth-order  functions  by  actual  functions.  This  yields 

S(x,t|r;x* ,t* |r' )  = 

r  .3  rT  .  i  ,  3G(x.t |r;y.3 is)  3?(x* ,t ' lr* ;y,s|s) 

Jd  y /  ds  U  (x,t|t;y,s|s)  - rl  ! - ^ — jr-*" - 

r  ,3  r1,  ,  i  ix  3G(&,t  |r;Y,s  |r)  3T(x*  ,t'  lr'  ^y.s  |r) 

tjd  yj  ds  U. .(5,t|t;y,s|s)  - ™  ! - ^ |  1 — 


+jd3yf  ds  U.  .(x,t  ir;y,s  |s)G(x' ,t*  |r';v,s|s)  — S ~ 
13  ~  *  '1 


♦Jd  yj  ds  1).  .(x,t  |t;y,s  !s)G(x‘ ,t'  |r’;y,sjr')  - a?x  1  dy 

-t  ~  ~  i  M 


(5.9) 


ri(x,t  i'rjx'  ,t'  jr' )  = 

3  r' 

Jd  yj  ds  U.j(x,t |t;y,s is) 


3G(&,t |r;y,s lr* )  3G(v,s |r* ;&* ,t*  jr* ) 

°X.  hYi 


U^d.tltjy.sU)  ^(^.tlr^.sU) 
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rig.  3  helps  in  understanding  (5.9)  and  (5.10).  The  first  two  terra* 
on  the  right-hand  side  of  (5.9)  arise  fro*  the  iteration  expansion  of  the 
#(x,t|r)  factor  in  (5.5).  The  first  tens  is  associated  with  the  evolution 
fro*  (tg.tg)  to  (r,r)  along  the  diagonal  in  Fig.  3a  and  the  second  term  with 
the  evolution  along  the  vertical  free  (r,r)  to  (t,r),  The  final  two  terms  in 
(5.9)  arise  fro*  the  ♦(x'.t'lr* )  factor  in  (5.5).  The  first  term  on  the  right- 
hand  side  of  (5.10)  is  associated  with  propagation  of  the  imposed  perturbation 
along  the  vertical  froai  (t*,r')  to  (r* ,r* )  in  Fig.  3b,  the  second  term  is 
associated  with  diagonal  propagation  fro*  (r*,r* )  to  (r,r),  and  the  third 
ter*  with  vertical  propagation  fro*  (r,r)  to  (t,r). 

For  further  clarification,  the  derivation  of  (5.19)  will  now  be  traced 

A 

in  aare  detail.  An  exact  expression  for  the  factor  G(x,t ir;x* ,t* |r' )  in 
(5.6)  is 


This  follows  i**ediately  fro*  the  equations  of  action  and  the  definition  of 
The  lowest-order  iteration  contribution  to  H  is  obtained  by  replacing 
u(x,t|t)  with  u^^x^tjt)  in  (5.6),  replacing  the  quantities  on  the  right- 
hand  side  of  (5.11)  with  zeroth-order  quantities,  inserting  this  result  in 
(5.6),  and  averaging.  Changinr  ail  zeroth-order  quantities  to  actual  quantities 
in  this  result  then  gives  (5.10). 
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Equations  ( 5. 1 )— ( 5.4) ,  (5.9),  (5.10),  (4.2),  and  (4.4)  are  a  complete 

set  which  determine  t(x,t ir;x' ,t ' | r’ )  and  G(x,t I r;x* ,t * | r' )  for  all  argument 

values  when  the  initial  function  f(x,tQ jt^ ;£* ,tQ jtQ )  is  specified  and  the 

Eulerian  velocity  covariance  LL.(x,t lt;x' ,t' S t • )  is  known  for  all  argument 

values.  Setting  t  =  r  ana  t *  =  r'  in  these  equations  gives  a  reduced  set, 

containing  only  Eulerian  quantities,  which  is  identical  with  the  Eulerian 

direct-interaction  equations  previously  obtained^  for  the  scalar  field. 

The  direct-interaction  equations  may  be  characterized  in  two  ways.  First, 

the  approximations  for  S  and  H  represent  the  summation,  to  all  orders,  of 

certain  well-defined  infinite  subclasses  of  terns  from  the  formally  exact 

iteration  expansions.  Second,  the  final  statistical  equations  are  exact 

statistical  equations  for  a  model  dynamical  system  which  has  some  important 

properties  in  cannon  with  the  actual  svsten.  These  matters  have  been  explored 

£ 

in  detail  in  the  pure  Eulerian  case. 

The  exact  statistical  quantities  satisfy 

jG(x,t !r;x* ,t! !r* )d3x  =  1  (5.12) 

and 

Js(x,t !r;x,t lr)d3x  =  0.  (5.13) 

For  t  =  r,  t *  =  r' ,  (5.12)  states  that  an  initially  introcucec  quantity  of 
scalar  fielc  substance  is  conserved  Dy  both  convection  and  molecular  diffusion. 
For  t  4  r,  t'  i  r',  it  states  that  t.ie  totai  Quantity  of  scalar  substance  is 
indepencent  of  whether  the  fiela  is  describee  in  Eulerian  or  Lagrangian 
cocrcinates.  (Cf.  (2.7).)  Eq.  (5.13)  expresses  conservation  of 

J[v(x,t  lt)]2d3x 

bv  the  convection  process  and  also  implies 


3  3 

/*(x,t ir;x,t ir)a  x  =  0, 


(5.14) 
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which  sta^-ss  that  the  aean-square  quantity  of  scalar  substance  is  independent 
of  whether  ascription  is  in  cuierian  or  Lagrangian  coordinates.  An  additional 
exact  property  is  that  when  *  -  0  the  functions  ?(x,t  jr;x'. ,t'  |r' )  and 
G(x,t |r;x' ,t* |r* )  are  independent  of  r  and  r'.  This  is  because  the  scalar 
concentration  is  constant  along  the  particle  paths. 

All  the  above  conservation  and  invariance  properties  survive  exactly 
in  the  direct -interact ion  equations.  Eos.  (5.12)  and  (5.13)  nav  be  verified 
by  integrating  (S.9)  and  (5.10)  over  space,  transferring  the  results  by 
partial  integration,  and  noting  that  the  U.  factors  are  solenoidal  in  i. 

The  independence  upon  r  and  r '  when  «  =  0  nay  be  verified  as  follows.  If 
all  the  quantities  on  the  right-hand  sides  of  (5.9)  ana  (5.10)  are  independent 
of  the  Measuring  tines  (the  tines  following  the  vertical  bars),  then  S  and  H 
are  independent  of  r  and  r*.  If  so,  and  «  *  0,  (5.1)  and  (5.2)  rive  identical 
Changes  of  the  functions  along  the  diagonal  of  the  (t,r)  plane  and  oarallel 
to  the  t  axis.  Therefore  the  change  parallel  to  the  r  axis  is  zero,  and  no 
dependence  on  Measuring  tine  is  generated  by  tne  equations.  Following  this 
argunent  in  further  detail  denonstrates  the  conplete  r  and  r'  independence. 

A  further  property  of  the  exact  statistical  equations  which  survives  in 
the  direct-interaction  approximation  is  the  existence  cf  equipartition  equili¬ 
brium  solutions  in  the  case  *  =  0.  If  f  is  appropriately  normalized,  these 
solutions  satisfy 

G(x,t ir;x' ,t* ir' )  =  G(x,t |t;x* ,t' It* ), 

f(x,t|r;x* ,t* |r*)  =  G(x,t |r;x* ,t* ir' ),  (5.15) 

when  r  >  r* .  The «qui  partition  solutions  can  be  inferred  for  the  exact 
equations  fron  the  existence  of  Liouville  and  fluctuation-dissipation  theorems 
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(see  Appendix).  No  condition  on  the  velocity  field  other  than  incompressibility 
and  vanishing  normal  component  on  the  boundaries  is  implied.  The  equinartition 
solution  is  physically  unrealizable,  but  it  is  important  because  it  expresses 
the  tendency  of  the  convection  to  produce  ever-sharper  gradients  of  tne 
scalar  field.  If  complete  equilibrium  were  ever  achieved,  than,  according  to 
(5.15)  and  (4.4), 

f(x,t |t ;x' ,t J t)  =  53(x-x'),  (5.16) 

which  indicates  zero  correlation  length  of  the  scalar  field. 

The  compatibility  of  (5.15)  with  the  direct -interaction  equations  when 
k  =  0  may  be  demonstrated  by  using  (5.15)  to  replace  all  ¥  functions  with 
G  functions  in  (5.1)-(5.4),  (5.9),  and  (5.10),  using  (4.2)  when  needed.  If 
now  a  partial  integration  over  y  is  performed  on  the  third  and  fourth  terms 
on  the  right-hand  side  cf  (5.9),  using  the  fact  that  the  factors  are 
solenoidal  in  j ,  then  it  is  easily  seen  that  cancellations  occur  in  such 
fashion  as  to  nake  (5.1)-(5.4)  identical  equations  for  G(x,t |t;x’ ,t' jt* ). 

All  the  above  consistency  properties  of  the  direct-interaction  equations 
can  be  predicted  frao  the  existence  of  a  rodel  representation  of  the  aoproxi- 
aation.  A  final  i no art ant  property  implied  by  the  existence  of  the  node! 
representation  is  that  the  covariance  f(x,t |r;x' ,t' |r’ )  predicted  by  the 
direct-interactioG  equations  is  realizable  as  the  average  over  some  possible 
ensemble  of  fields  and  therefore  satisfies  all  the  realizability  inequalities 
which  a  covariance  must. 


\ 


4 
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6.  LHDI  EQUATIONS:  SCALAR  FIELD 


Suppose  for  now  that  the  spatial  domain  is  infinite  and  the  fields  are 
statistically  homogeneous.  Consider  the  function 
S(x,t(t;x',t|t)  = 

Jd3y/  ds  Uj  .(x#t  | t;y ,s j  s)  G(x,t|  t;y,s|  s)  5 } 

~  i  1  ~~  yi 


♦  G(x'.t|t;y,s|s) 


(6.1) 


given  by  (S.9).  The  Fourier  transform  of  (6.1)  with  respect  to  x-x'  yields 
the  direct-interaction  approximation  for  the  transfer  of  mean-square  scalar 
substance  between  the  wavevector  modes  of  the  scalar  field.  The  integration 
over  s  in  (6.1)  expresses  relaxation  and  memory  effects  in  the  turbulent 
flow.  Clearly  the  magnitude  of  the  integrals  is  dependent  on  the  Eulerian 
correlation  times  of  the  velocity  and  scalar  fields.  These  correlation  tiroes 
play  the  role  of  effective  tiroes  for  renemfcering  contributions  to  the  transfer 
function  from  dynamical  interactions  during  the  oast  history  of  the  fluid. 

The  appearance  of  the  Eulerian  correlation  tiroes  in  (6.1)  gives  rise  to 
a  serious  flaw  of  the  direct-interaction  aporoximation  in  the  description  of 
convection  effects.  The  trouble,  which  is  discussed  in  detail  in  an  earlier 
paper, *  can  be  stated  as  a  failure  of  the  direct-interaction  approximation  to 
preserve  the  invariance  properties  of  the  exact  equations  under  a  random 
Galilean  transformation.  Supoose  that  u(x,t)  is  augmented  by  an  addition  v 
which  is  constant  in  space  and  time,  statistically  indooendent  of  u(x,t)  at 
any  instant,  and  Gaussianly  distributed.  This  means  that  the  systems  in  the 
ensemble  are  subjected  to  uniforms  translations  that  differ  randomly  from 
system  to  system.  The  translations  do  not  affect  the  rate  of  transfer  of 
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scalar  substance  among  the  wavevector  modes  in  any  of  the  systems  of  the 
ensemble,  and  therefore  it  is  clear  that  S( x,t | t ;x* ,t | t )  must  be  invariant 
under  the  random  Galilean  transform, ation.  But  (6.1)  is  not  invariant,  because 
the  i-ulerian  correlation  times  are  net  invariant.  The  latter  deoend  on  how 
fast  the  velocity  y  sweeps  the  fluctuations  in  the  v  fiexd  past  fixed  ob¬ 
servation  ocints  in  soace. 

The  following  question  nay  be  posed:  Is  it  possible  to  alter  the 
direct-interaction  equations  so  as  to  incorporate  invariance  under  randc*® 
Galilean  trams format ion,  without  giving  uo  any  of  the  conservation,  invariance, 
and  equilibrium  properties  exhibited  in  Sec.  5?  Investigation  has  so  far 
indicated  that  the  objective  can  be  achieved  only  in  part.  There  appears  to 
be  a  unique  prescription  for  altering  the  direct -interact ion  formulas  for  S 
and  rl  so  that:  a)  The  conservation,  invariance,  and  equilibrium  properties 
of  Sec.  5  are  all  preserved,  b)  Exact  invariance  to  random  Galilean  trans¬ 
formation  is  realized  in  the  case  <  =  0.  c)  'When  the  velocity  field  con¬ 
sists  wholly  of  a  Gaussianly  distributed  uniform  velocity  v,  the  results  for 
f(x,t jr;x' ,t ' | r* )  and  G(x,t |r;x' ,t' |r' )  are  exactly  correct,  without  restriction 
on  <.  The  prescription  is: 

(1)  Write  U. .(x,t jt ;y ,s Is)  in  (b.9)  and  (5.10)  as  U. .(x,t lt;v,s is)  so 

.i]  ~  ~  i]  '  - 

as  to  exhibit  explicitly  the  fact  that  it  is  solenoidal  in  j. 

(2)  Then,  everywhere  on  the  right-hand  sides  of  (5.9)  and  (5.10), 
change  each  labeling  time  s  to  a  t.  (The  labeling  times  are  the  tine 
arguments  which  precede  the  vertical  bars.) 
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£q.  (6.1)  is  now  replaced  by 


,t| t)  = 


The  s  integration  in  (6.2)  is  over  Lagrangian  histories*  and  it  is  now  the 
correlation  tines  which  detemine  ivow  the  transfer  function 
re  ue  labors  p*st  dynanical  events.  More  generally*  aoplication  of  the  pre¬ 
scription  to  (5.9)  and  (5.10)  gives 


S(*,t|r;x' ,t* |rf)  = 
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In  the  second  term  on  the  right-hano  side  of  (6.3)  and  the  third  term  on  the 

right-ha;. d  side  of  (6.4),  the  y  integration  ;  as  been  performed  bv  using  (4.4) 

and  the  fact  that  the  factor  U ,^( x,t I t ;y ,t ! s)  is  solenoidal  in  i.  Eos.  (6.3) 

i]  “  - 

and  (6.4)  will  be  called  the  LriDI  ( Lagrangian-histcrv  direct-interaction) 
approximation  hereafter. 

The  invariance  of  the  LnEI  equations  to  random  Galilean  transformations 
will  be  discussed  in  a  moment.  First,  it  is  aasy  to  see  that  the  conservation, 
invariance,  and  equilibrium  properties  noted  in  Sec.  5  all  survive.  For  the 
unaltered  direct -interact ion  equations,  tnese  properties  depend  upon  tne 
cancellation  or  identity  of  particular  integrals  in  the  expressions  for  S  and 
n.  The  alteration  does  not  change  the  limits  of  the  integrals,  and  identical 
integrands  are  subjected  to  identical  changes.  Thus  the  demonstration  of  all 
the  properties  goes  through  for  the  LHOI  equations  just  as  before.  It  is  im¬ 
portant  to  note  that  if  rule  (1)  of  the  alteration  prescription  were  not  im¬ 
posed  tne  equilibrium  property  would  not  survive.  This  property  reauires  that 
the  velocity  covariances  which  appear  be  solenoidal  in  j ,  and  U„(x,t  !t;y,t  Is) 
need  not  be  solenoidal  in  j  for  s  /  t. 

The  Galilean  transformation  properties  are  conveniently  investigated  by 
writing  the  equations  in  tne  wave-vector  domain,  using  (4.7).  Eqs.  (5.1)-(3.4), 
(6.3),  anc.  (6.4)  are  replaced  by 


(3/dt 

♦  <k2)'f(k;t  it  ;t’ 

Ir') 

=  S(k  ;t 

1 1  ;t ' 

ir'). 

(6.5) 

0/3t 

♦  kk*  )G(k;t 1 1 ; t  * 

|r* ) 

=  ii(k;t 

It;t* 

Ir'), 

(6.6) 

3'f(k  ;t  ir;t’  jr' 

)/3t 

=  3(k;t 

i  r;t ' 

ir’). 

(6.7) 

3G(k ;t ir;t ’  |r’ 

)/3t 

=  rt(k;t 

ir;t  * 

ir’). 

(6.3) 

and 


I 


-2**  - 


S(k;t jr;t* jr* )  * 


rA  -  S 

*kik3T>.q^  dS  ^^ItStUXKpstlritlsmkjtVitls) 

*  ”  to 

-k  k  »(k;f|r';t!r)Jd3,/  ds  U.?<d;t I t;t is) 

J  _  1  I 


*  Q 

<**  |s)G(ik;t*  lr'jtis)t(p;tlr;t|s) 


:A  f1'  ..  S 


+kikj^P.q^  Uij(a;tlt;t!s)G(ii;t’lr,;tlr’)np;t|r;tlrM,  (6.9) 


-  r 


H(k;t  r;t' |r')  = 


‘kikj^P,a/t,ds  Uij(i;t  lt;t  is)G(p;t  Ir;t  lr’  )G(k;t  |r*  ;t • jr' ) 


^  i  S 

“kikjlp^q/  Ui.(q;t|t;t|s)G(p;tir;t!s)G(k;t|s;t’ lr’) 


_  t 

-k.k  G(k;t|r;t’ |r’)/d  qj  ds  U.T(q;t |t;t !s). 


In  these  equations  the  operator  notation 


Ip#q  5  /d3pd3q63(k-p-q) 


(6.10) 


(6.11) 


and  the  solenoidal  properties  of  an*  used. 

Suppose  that  the  random  uniform  velocity  v  is  statistically  isotropic. 
(This  is  an  inessential  restriction.)  Then  the  exact  effects  of  the  random 
translation  on  G  and  f  are  expressed  by 

[G(k;t|r;f  ir’)Jv  =  exo[-  y  v2k2(t-t’  )?][G(k;t  |r;f  jr’  ) ] 

Cnk;t|r;f  |r’)]y  =  exp[-  y  v2k2l t-t * )2]lf(k;t lr;t’  |rf  )j  , 


(6.12) 


f 
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where  v  is  the  res  value  of  any  vector  component  of  v,  [  ]  denotes  a  value 

J  **  0 

with  v  absent,  and  [  ]  denotes  a  value  with  v  present.  Eq.  (6.1?)  follows 
frees  the  Gaussian  distribution  of  v^  and  the  fact  that,  according  to  (2.9)- 
(2.11),  the  scalar  Fourier  coefficients  in  the  individual  systems  of  the 
ensemble  transform  according  to 


[f(k,t|r)]v  =  exp[-iv*k(t-t0)][it»(ik,t  !r)]Q. 

(Cf.  tef.  1.) 


(6.13) 


*,  .  (6.5)-(6.3)  and  (6.12)  imply  that  S  and  H  transform  according  to 
[S(_;t|r;f  |r‘)]v  =  expl-  A  v2k2(t-t’  )2]{[S(k;t  |r;f  |r* )  }Q 

-  v2k2(t-t*  )[t(lc;t  |  r;t '  |r*  )]Q), 

[rf(kjt jr;t' |r' )]y  =  exo[-  A  v2k2(t-t ’  )2]{  [Hi k ;t  | r;t *  j r*  ) ]Q 

-  v2k2(t-t')CG(k;tir;t’ |r')30).  (6.14) 


The  effect  of  the  transformation  upon  the  velocity  covariance  itself  is 


tU„(q;t  |r;t’  |r’ )  jy 


6ijfiJ(q)vQ  ♦  exp[-  y  v2k2(t-t’ )^3[U.t(q;t ir;t' |r’ )]Q. 

(6.15) 


Now  take  <  -  0.  Then  the  LriDI  equations  yield  *( x,t I r;x* ,t * I r* )  and 
G(x,t |r;x' ,t' |r' )  values  which  are  independent  of  r  and  r' ,  in  agreement  with 
the  exact  functions.  Using  this  fact  (noted  above),  it  is  verifiable  by  direct 
evaluation  that  the  substitution  of  (6.12)  and  (6.15)  into  (6.9)  and  (6.10) 
gives  (6.14).  This  demonstrates  the  desired  Galilean  invariance  prooerty. 

It  is  important  to  remark  that  this  result  would  not  be  obtained  if  any  or 
all  of  the  labeling  times  s  were  changed  to  any  value  other  than  t,  in  forming 
the  LiiDI  equations. 
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Next  take  <c  >  0  and  let  the  velocity  field  consist  solely  of  a  random 
uni for®  field: 

Uij^l*tlr*t'  ir' )  =  6ij63(3.)vO*  (6.16) 

Direct  solution  of  the  equations  of  motion  in  Secs.  2  and  3  then  yields  the 
exact  values 

T(k;t  |r;t'  |r' )  *  exp[-  y  v2k2(t-t * )2  -  Kk2(r+r*  )3*(k;tQ ! tQ  jtQ  ltQ) , 

G(k;t |r;t*  fr' )  =  expt-  j  v2k2(t-t ' )2  -  *k2(r-r')]  (r  >  r*).  (6.17) 

Direct  evaluation  shows  that  the  values  (6.16)  and  (6.17)  are  also  exact 

solutions  of  the  LhDI  equations  (f  ,)-(6.10). 

In  the  case  of  a  general  isotrooic,  homogeneous,  and  statistically 

stationary  velocity  field,  a  Peclet  number  B  may  be  defined  by 

B  =  (6.18) 

where  vQ  is  now  the  total  rms  value  of  any  vector  component  of  velocity  and 

t  is  an  appropriate  correlation  scale- length  for  the  velocity  field,  say  the 
7 

integral  scale.  In  the  cases  of  uniform  translation  treated  above,  the 
correlation  scale  is  infinite  and  so,  therefore  is  b.  When  3  <<  1,  the 
turbulent  convection  represents  a  small  perturbation  on  the  molecular  diffusion 
process,  and  the  Lulerian  and  Lagrangian  histories  of  the  scalar  field  differ 
inappreciably.  In  this  case  the  unaltered  direct -interact ion  and  LHDI 
aoproxiaations  give  nearly  the  same  results.  A  formal  statement  of  this  fact 
can  be  made  by  noting  that  the  iteration  expansion  discussed  in  Sec.  5  is 

y 

actually  an  expansion  in  powers  of  b.  The  exact  functions  S  and  H  are  tnen 
power-series  in  B,  with  the  lowest-order  terms  linear  in  3.  If  the  tiireev- 
interaction  and  LHDI  values  for  S  and  H  are  expanded  in  powers  of  B,  the  terms 
linear  in  b  agree  with  the  exact  expansion.  In  the  higher  powers,  the  two 
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approximations  differ  from  each  other  and  from  the  exact  exoansion.  In  the 
asymptotic  situation  of  convection  by  a  random  uniform  velocity,  treated 
above,  the  LnDI  values  agree  with  the  exact  expansion  to  all  orders,  while 
the  unaltered  direct-interaction  values  give  the  higher  orders  with  rapidly 
deteriorating  accuracy.  (See  Ref.  8) 

It  should  be  stressed  that  one  consistency  property  of  the  direct- 
interaction  equations  has  not  been  shown  to  survive  in  the  LHDI  equations. 

No  model  representation  sc  far  has  been  found  for  the  latter,  and  consequently 
it  cannot  be  asserted  that  the  function  *(x,t Jr;x’ ,t5 |r* )  they  predict  is 
a  realizable  covariance.  This  matter  will  be  discussed  in  the  next  Section. 


7.  EXAMPLES  OF  THE  SCALAR  EQUATIONS 

Sec.  6  contained  a  demonstration  that  the  LHDI  equations  yield  exactly 
correct  V  and  G  values  for  convection  by  a  Gaussianly  distributed  uniform 
and  constant  velocity  field.  This  is  a  nontrivial  success.  The  unaltered 
direct -interact ion  equations  predict  a  peculiar  wave-front  behavior  for 
G(x,t |t;x' ,t' j t * )  in  the  same  situation,  while  cumulant -discard  approximations 
or  truncations  of  the  iteration  expansion  lead  to  disaster.-'*  The  success 
here  gives  a  degree  of  reassurance  concerning  the  absence  of  j.  demonstrated 
model  representation  for  the  LHDI  equations.  The  question  of  model  representa¬ 
tion  has  assumed  importance  because  the  curulant=discard  aDproxim«itions, 
truncations  of  the  perturbation  expansion,  and  certain  higher  infinite  partial 
sunnatio^s  of  the  perturoation  expansion,  all  of  which  have  no  model 
representations,  have  been  found  to  lead  to  grossly  unphysical  behavior. 

9 

This  has  taken  the  form  of  negative  wave-vector  spectra  in  the  scalar  problem 
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and  in  Xavier-Stokes  turbulence  dynaoics.*0  With  each  of  these  approxi nations , 
the  unphysical  behavior  is  apparent  in  the  problem  of  convection  by  a  random 
uni fore  velocity  field.6*8 

It  has  been  noted  that  G  and  f  are  independent  of  the  aeasuring  tines 
when  k  s  0.  In  this  case  the  G  equation  given  by  (5.2)  and  (6.4)  reduces  to 


G(*#t»*' •*')  *  377  >  /  u.^(x,t|t;x,t js)ds  , 


3x.  [“ 


G(x,t' ,t* )  =  6  (x-x* ), 


(7.1) 


where  G(x,t;x‘«t*)  is  written  for  G(x,t |t;x* ,t' |t' )  and  (4.4)  has  been  used. 


If  IL  .  has  the  isotropic  homogeneous  forte 


Ui.(x,t|r;x,t|s)  =  6_U(t|r;t  s). 


(7.2) 


then  (7.1)  has  the  immediate  solution 


G(x#t;x\t‘)  =  [2«o1(ttf>3"3/2exp[-|x-x»|2/2ol(t#t*)3, 

where  t  s 

a  (t,t*)  =  2 /  ds/  dr  U(s|s;s|r). 

t*  t* 


Eq.  (7.4)  yields 

(U-x*  |2>  s 


(7.3) 

(7.4) 


(7.5) 


These  results  nay  be  compared  with  exact  formulas  obtained  by  the  method 
11  5 

of  Taylor,  *  The  displacement  of  a  particle  from  x*  to  x  in  the  interval 
from  t*  to  t  is,  by  definition. 


hence. 


where 


x-x*  =  J  u(x,tjs)ds. 

t'~~ 

(|x-x* i2>  =  30(t,t* ), 


o(t,t*)  =  /  dr/  ds  U(t!s;t|r). 
t*  t* 


(7.6) 


(7.7) 
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If  the  statistical  dependence  among  the  values  of  u(x,tis)  has  an  effectively 
finite  range  in  t-s,  then  a  central- limit  argument  inplies  that  for  very  large 
t-t'  the  function  G  is  Gaussian;  that  G  is  given  by  (7.3)  with  o  reolaced 
hv  o.  On  the  other  hand,  Roberts5  has  found  that  for  very  snail  t-t’  the 
function  G  is  again  Gaussian.  Supoose  now  that  the  turbulence  is  statistically 
stationary.  Since  u(x,tjs)  and  u(x,tir)  are  measured  along  the  same  particle 
path,  it  follows  from  the  stationarity  and  honsogeneity  that  U(tjs;t|r)  can 
depend  only  on  js-rj.  Then  U(t|s;t|r)  =  U(sls;s|r)  and  c(t,t’)  =  o^t.t’). 

In  this  case  the  LnDI  result  for  G(x,t;x’,t')  is  asymptotically  exact  for 
very  o-iall  and  very  large  t-t',  while  the  result  for  (lx-x'l^>  is  exact  for 
all  t-t'. 

The  spectral  transfer  equation  yielded  by  (6.S)  and  (6.9)  is 


O/at  ♦  2«k  j¥(k ;t  jt  ;t  j  t )  =  2k  k  j  J  U.*J( q;t !  t  ;t  I  s )ds 

**  1JP  *9 


*  [G(k;t|t;tjs)V(p;t|t;t|s)  -  G(  p;t  1 1  ;t  i  s )  ’»'( k  ;t  i  t  ;t  |  s)  ].  (7.8) 

When  <  =  0,  the  G  and  t  functions  are  independent  of  iabelinr  tines,  and  (7.8) 
takes  the  sicole  form 


A  r 

3T(k;t)/at  s  2k.k4J?  [»(?it)  -  f(k;t)]/  u.^(q;tit;t|s)as. 


(  7.9) 


where  *(k;t)  =  *(k;t  |t;t  |t)  and  (4.4)  is  used.  The  x-soace  equivalent  of 
( 7.9 ) .obtained  from  (6.2)  specialized  to  the  homogeneous  case,  is 


at 


3  ^(x-x'.t)  =  3  "l"^'  /  3.  .(x-x*  ;t|t;t  |s)ds, 

x.  X.  tQ  13  ' 


(7.10) 


where 


3.  .(x-x*  ;t  It  ;t  |s)  =  2fJ .  S.(  x,t  j  t  ;x,t !  s )  -  U  5(x,t  it;x*  ,t  is)]  (7.11! 


and  'f(x-x’.t)  =  Y(x,t  J  t;x’  ,t ,  t ).  When  the  turbulence  is  statistically 
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stat ionary  as  well  as  homogeneous. 

U  ^(x,t|t;x,t|t)  =  U.^x*  ,t|s;x*  ,t|s), 

A  J  A  j 

in  accord  with  a  remark  Bade  above.  In  this  case.  B. .  is  the  variance  tensor 

il 

of  the  velocity  difference  u(x,t|t)  -  u(x* ,t|s). 


8.  LHDI  EQUATIONS :  VELOCITY  FIELD 

The  direct-interaction  approximation  and  Lagrangiar-history  alteration 
go  through  for  the  velocity  field  in  close  analogy  to  "  ;  scalar  case.  The 
analysis  in  this  Section  will  be  restricted  to  reflection- invariant  homo¬ 
geneous  turbulence.  More  general  flows,  with  nonzero  Bean  fields,  can  be 
treated  with  the  netnods  of  Kef.  3.  Assume  that  at  tine  tQ  the  homogeneous 
Eulerian  field  is  Gaussianly  distributed  and  divergence less.  Then  the 
different  Fourier  components  are  statistically  independent  and  the  initial 

s 

distribution  is  fully  specified  by  the  function  U^.(k;tQ jtQ;t0 jt^).  The 
exact  statistical  equations  are 


(a /at 

*  vk2)U.  (k;tjt;c* 
m  - 

ir') 

=  S.  (k;t  jt ;t ' 

k’>. 

(8.1) 

(J/at 

♦  vk2)G.  (kjtjtjt' 
in  -  *  * 

|r») 

=  H.  (k;t|t;t* 
m  -* 

k’>. 

(8.2) 

JU.  (k;t jr;t*  Jr* 
in  - 

)/3t 

=  «.  (k;t  jr;t  * 
in  - 

k’>. 

(8.3) 

3G .  (k;t |r"t' |r' 

)/3t 

-  it’.  (k;t  jr;t* 
in  - 

k'). 

(8.4) 

where  S.  ,  ri.  ,  M.  ,  and  N.  are  the  Fourier  transforms,  with  respect  to  x-x’ , 
in*  in*  in*  m 


SiQ(x,t  i t ,t*  j r* )  = 


-  i  Ku^x.tjtJu^Cx.tltJu  (x\t*|r*)>,  (8.5) 

l  X  1  ^  ■  o 


of 


-31- 


i  ^ 

tv  t  1 1  •  v  *  r  *  i«*  \  —  /  c  \  /  ..“  /  w  *-  i  «>  i  x*  t ..  ^  i  ^  _■  .  i  i  •  A  /  *  ,  » 

in'*’-  ’=  ••  ’•  '  w^jnx~’u’  k's  f"  ir//* 

Kln(x,tir;x',t'lr')  =  -(u^x,tlt)Un<x\tMr')3u.(x,t;r)/3x]i>,  (6.7) 

Hin(x,t|rix'.f!r-)  =  -^(*,t|t)JC.n(x.tlr;x',t1lr')/3xi[> 

.<aSx,t!t;x,,t‘irM3u;(x,t!r)'3x  ).  (e.6) 

The  direct-interaction  approximations  for  S.  .  n.  *.  .  and  N  are 

m’  ms  m*  'in 

formed  by  the  same  rules  as  in  the  scalar  case:  Iteration  expansions  for 
these  quantities  are  constructed  as  functional  power  series  in  the  zeroth— 
order  covariances  and  Green’s  functions.  The  latter  satisfy  the  zercth- 
order  equations  formed  by  discarding  the  right-hanc  sices  of  (3.18),  (3.21), 
(3.24),  and  (3.25).  Then  the  leaning  terms  only  are  retained  in  the  expansions, 
and,  in  them,  all  zeroth -order  quantities  are  repxacad  bv  actual  quantities. 

The  prescription  for  obtaining  the  LnDI  approximations  tor  the  velocity 
field  is  simply: 

On  the  right-hand  sides  of  the  direct-interaction  expressions  for  S  . 

in’ 

riin»  Mia»  3041  «in  change  every  labeling  time  s  to  a  t. 

A  rule  like  rule  (1)  of  Sec.  6  is  not  needed  here  because  the  recuired  S 

superscripts  are  autonati rally  supplied  as  a  result  of  writing  the  equations 

of  motion  in  the  forms  (3., 18)  and  (3.21). 

The  final  results  for  the  L'-tDI  approximations  to  S.  ,  ri  ,  ?.  .  ana  »'  are 

m*  in*  m*  'in 

Sin(~ltit;t,lr,)  = 

"Pijc(^)^,q?bca(f  }f  ^•b(p;tit;tis)U^(q;t!t;t,s)Un^(k;t’  <r’;t|s) 

t0 

r  * 

+  2  Pijm(-)?abc(-)^,qjr  45  Gna(-;t'  I1"’ ^  ls>u^(?;t  !t;t  |S)U^(q;t  jt;t  |s) 

"  t0 

A  t' 

^ijm^^p.q rafr,Gr’  Gnb(^;t’  ;r* ;t  |r’  )Jjt(L:t  !t;t  !r’  )jra(  l;t  |t;t  |s)*  (3’S) 


-33- 


where 

P.  .  (k)  =  k  P.  .(k)  ♦  k.P.  (k).  (8. 13) 

1 1  it.  i]  -  i  id  ~ 

The  paths  of  evolution  and  propagation  involved  in  these  equations  are 
indicated  in  rip.  4.  The  solenoidal  properties  ana  (4.9)  have  been  used  in 
writing  some  of  the  terms. 


Lqs.  (8.1)-(8.4),  ( 8.9 )-( 8. 12) ,  and  (^.9)  form  a  complete  set  which 

determine  the  full  functions  U.  (k;t jr;t* |r* )  and  G.  (k;t ir;t* |r' )  from  the 

in  in  'w 

prescribed  initial  values  U.  (k  ;t.  |t_ ;t, Jt_).  As  in  the  scalar  case,  thev 

in  ~  j  u  u  u 

preserve  important  oroperties  of  the  exact  equations  of  motion:  Conservation 
oi  energy  by  the  nonlinear  interaction;  tne  invariance  property 


3 

3t 


Jo. .(x,t ir;x,t ir)dJx  =  3; 

i] 


(8.14) 


ana  the  existence  of  formal  equipartition  equilibrium  solutions  of  the  form 


U.^(k;t  lr;t*  |r* )  =  3_(k  ;t  jr;t  ’  lr’ )  (r  ±  r* )  (8.15) 

in  the  case  v  =  0.  A  property  of  the  aprroxir.at ion  which  has  no  analog  in 

the  scalar  case  is  that  the  incompressibility  of  the  bulerian  velocity  field 

C  i 

is  preservea:  ( k ;t I t ;t 1 | r ’ )  vanishes  tor  all  argument  values  if  it  is 
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ser©  for  t  =  t*  =  r'  =  tQ. 

All  of  these  consistency  properties  can  be  inferred  for  the  unaltered 
direct-interaction  equations  fro*  the  existence  of  a  model  representation. 
As  in  the  scalar  case,  they  can  be  verified  for  the  LHDI  approximation  by 
direct  examination  of  the  final  statistical  equations.  The  cancellations 
of  terms  in  (8.9)-(8.12)  associated  with  the  conservation,  invariance,  and 
equipartition  properties  can  be  demonstrated  by  using  the  solenoidal 
property  with  respect  to  indices  beneath  S  suoerscripts  and  noting  the 
geometrical  relations 


k  =  P  ♦  3  #  P..(k)  =  P.  .  ( p )  ♦  P  ..(q),  (8. 16) 

—  --  _  13m  -  jmi  £  mi 3  2. 

The  preservation  of  incompressibility  can  be  verified  by  tracing  the  oath 

of  evolution  of  U^(k;t  jt;t*  Jr’ )  and  IL  .(k;t jr;t' jtJ )  and  noting  that  the 

presence  of  projection  operators  P  or  superscripts  S  in  (8.9)  and  (8.11) 

precludes  the  generation  of  any  increments  which  are  nonsolenoidal  in  i. 

Eq.  (4.2)  is  needed  in  the  demonstration.  [Some  of  the  S  superscripts  in 

(8.9)-(8.12)  are  superfluous  because  of  (4.10).  They  are  retained  in  order 

to  make  it  easier  to  trace  the  way  in  which  the  various  terms  arise  from  the 

perturbation  expansion.] 

Equations  ^i.^l)  and  (3.25)  are  not  fully  invariant  to  Galilean  trans¬ 
formation  as  they  are  written.  To  obtain  invariance,  the  respective  terms 

(8.17) 

must  be  added  to  the  right-hand  sides,  where  v  is  the  uniform  translation 
velocity.  This  is  a  formal  device  which  insures  that  "impossible”  perturba¬ 
tions  (see  Sec.  3)  are  translated  along  with  everything  else.  The  added 
terms  make  no  change  in  the  evolution  of  u(x«tjr)  when  the  initial  .Suierian 


c  c 

•Wu.  (x,t  it' ,  -W  G.  (x,t  it  ;x*  ,t' ir* ) 
~  1  1  -  ~x  m  -  ~  ’  1 
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field  is  divergence less.  With  this  change,  it  is  easily  verified  that  the 
exact  statistical  functions  transform  as  follows  under  the  isotropic  random 
Galilean  transformation  considered  in  Sec.  6: 

[U.n(k;t|r;t*  |r*)Jv  =  exp[  -  ±  v2k2(t-f  )2J[U.n(k;t !r;f  |r» )  JQ,  (8.18) 

CG^n(^*t i r;t* lr* ) 3y  =  exp[  -  J  v2k2(t-t» )2)[Gin(k;t|r;t* Jr* ) ]Q.  (8.19) 

The  notation  here  is  the  same  as  in  Sec.  6.  A  further  transformation  law  is 

^Sin^*f  *qit;t'  >^v  =  *Xp^  ”  7  vok2(t”t*  )23CSin(k;p,q;tit*  )]q  *  °)» 


where  S^n(k;p,q;t ;t' )  is  defined  by 


(8.20) 


(8.21) 


Eq.  (8.20),  and  similar  transformation  laws  for  other  functions,  actually  is 
implied  by  (8.18). 

By  analogy  with  Sec.  6  it  might  be  expected  that  the  velocity-field 

L'.iDI  equations  exhibit  invariiince  to  random  Galilean  transformation  whenever 

v  =  0.  This  is  not  the  case.  Verification  of  Galilean  invariance  in  Sec.  6 

depended  on  the  fact  that  the  scalar  statistical  functions  are  independent  of 

measuring  time  when  «  =0.  However,  variation  of  the  velocity  along  the 

particle  paths  arises  from  pressure  as  well  as  viscous  forces.  Consequently 

Uin(x,t !r;x' ,t' jr’ )  can  vary  with  r  and  r*  when  v  =  0.  The  velocity-field 

LKDI  equations  exhibit  exact  Galilean  invariance  only  under  very  restricted 

circumstances:  when  v  =  0  and  the  spatially  nonuniform  part  of  the  velocity 

field  is  infinitesimal  in  comparison  to  v  .  In  this  case,  for  elapsed  times 

not  indefinitely  large  compared  to  (vQk)  \  variation  of  lKn(x»t lr;x' ,t' lr* ) 

and  G. .(x,t |r;x’ ,t’ Jr* )  with  r  and  r*  can  be  neglected  in  computing  S..,  etc. 
i]  —  ~  *j 

Then  it  can  be  verified  by  direct  substitution  that  the  LHDI  equations  yield 
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(8.18)-(8.20).  In  carrying  out  the  demonstration,  care  must  be  taken  tc 

include  in  the  LHDI  equations  all  the  extra  terms  induced  by  the  additions 

(8.17)  to  the  underlying  equations  of  motion.  A  ooint  of  consistencv  i< 

that  these  extra  terms  turn  out  not  to  affect  the  evolution  of  U  (k ;t |r;t ' ir' ) 

in  ~  ' 

°r  °f  Gin^*t®r*t*  !r*>»  provided  that  :j£n(i.it0!t0it0't0)  vanishes. 

The  invariance  Droperties  of  the  LHDI  equations  u* cer  random  Galilean 
transformation  may  be  compared  with  the  behavior  of  the  unaltered  direct- 
interaction  equations.  For  t  =  t*,  (8.20)  expresses  the  fact  that  tne  random 
uniform  velocity  does  not  affect  the  simultaneous  triple  correlations  induced 
by  interaction  of  the  various  nonzero  wave-vector  components.  In  contrast,^ 
the  unaltered  direct-interaction  approximation  predicts  a  spurious  aecav  of 
these  correlations  with  a  characteristic  tine  Mv  k)"1. 


9.  RLLATIOH  TO  TriE  KOLMOGOROV  THEORY 

When  tne  turbulence  is  isciropic,  U.  and  G.  take  the  forms 

in  1 « 


m 


U.  (k;t|r;t*  |r* )  =  P.  (k)U°(k;t!r;t* Ir* )  *  F.  (k  )'JL(k  ;t  I  r;t  *  i  r* ) , 
m  —  '  m  -  m  —  ’  ’  * 


m 


m 


m 


G.  (k;t  |r,t* |r* )  =  P.  ( k )GS( k ;t ! r;t ’ I r* )  ♦  F .  (k )Gu(k ;t lr;t * Ir' ) ,  (9.1) 

m  in-  in—  ’  ’  * 


m 


in 


and  the  LnDI  equations  of  Sec.  8  can  be  reduced  to  scalar  form.  In  ~>articular, 
(8.9)  gives  the  energy-transfer  function 

T(k,t)  *  **f  2k3 //^pqdodqj  ds[a( k  ,p  ,q)G~(k  ;t  1 1  ;t  i  s  )!JS(p ;t  1 1  ;t  i  s )U^(  q  ;t  1 1  ;t  ‘  s  ) 

t„ 


-Mk  ,p,q)G^(  p;t  lt;t|s)U°(q;t!t;t|s)US(k;tit;t!s)J, 

where  T(k,t)  is  defined  by 

Unk^s.  (k;t|t;t|t)  =  P,  (k)T(k,t). 
in  in 


(9.2) 


(9.3) 
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The  integration  J  J  in  19.2)  is  over  all  wavenumbers  p  and  q  which  can  form 

a  triangle  with  k,  and  the  geometrical  coefficients  a  and  b  are  the  same  as 

2 

in  the  unaltered  direct-interaction  formula. 

The  integration  over  the  past  in  (9.2)  involves  correlations  measured 
along  the  paths  of  fluid  elements,  in  contrast  to  the  unaltered  direct- 
interaction  app  oximation  where  the  time  correlations  are  taken  at  fixed 
points  in  space  (cf.  Sec.  6).  For  a  fluid  obeying  the  Navier-Stokes  equation, 
the  change  of  velocity  along  particle  paths  arises  solely  from  the  viscous 
and  pressure  forces.  If  the  terms  representing  these  forces  are  omitted  from 
the  equations  of  motion,  then  U  (k;t|t;tjs)  and  G  (k;tlt;t|s)  are  independent 
cf  s.  This  property  is  preserved  in  the  LHDI  equations:  If  the  viscous  and 
pressure  terms  are  dropped,  the  equations  for  propagation  of  the  statistical 
functions  parallel  to  the  t  axis  and  along  the  diagonal  in  the  (t,r)  plane 
become  identical,  and  there  is  no  variation  with  measuring  time  r.  Thus  (9.2) 
states  that  the  effective  relaxation  time  for  triple  correlations  of  simul¬ 
taneous  amplitudes  of  wave-vector  triads  is  determined  by  memory  and  decay 
times  associated  with  the  viscous  and  pressure  forces  encountered  along  the 
particle  paths. 

Now  consider  statistically  stationary  high-Reynolds-number  flow.  Suppose 
that  an  inertial  range  exists.  By  this  is  meant  a  range  of  wavenumbers  below 
which  lies  most  of  the  kinetic  energy  and  above  which  lies  most  of  the  mean- 
square  vorticity,  or  dissipation.  If  the  energy-spectrum  follows  a  newer  law, 
these  suppositions  require 

£(k )  «  k-n  (1  <  n  <  3) 

where  E(k)  =  2*k2U*"(k  ;t !  t  ;t  1 1 ) . 


(9.U) 
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Both  the  viscous  and  pressure  terras  in  the  Navier-Stokes  equation  involve 

velocity  gradients.  With  a  spectrum  of  the  forra  (9.4),  the  energv-range 

velocity  excitation  makes  a  negligible  contribution  to  the  mean-square 

acceleration  of  fluid  elements,  provided  that  the  flatness  factors  of  the 

statistical  distribution  of  the  energy-range  wavenumbers  are  not  too  large. 

This  implies  that,  if  k  is  in  the  inertial  range,  the  correlation  time 

t 

T k  s  J  U(k;t|t;t|s)ds/U(k;t|t;t|t)  (9.5) 

is  negligibly  affected  by  the  energy-range  excitation  and  is,  in  fact,  the 
order  of  the  reciprocal  rms  vorticity  associated  with  wavenumbers  of  order  k. 
That  is, 

Tk  %  (vkk)’1'  vk  =  CKE(k) J1/2.  (9.6) 

This  argument  carries  over  to  the  LHDI  equations  because  the  gradient  operators 
in  the  Navier-Stokes  equation  are,  of  course,  carried  through  the  perturbation 
analysis  that  underlies  the  approximation.  Tne  vorticity  associated  with 
wavenumbers  much  higher  than  k  does  not  contribute  appreciably  to  because 

_  3 

of  averaging  effects  over  spatial  volumes  of  order  k  and  times  the  order  of 

V 

The  above  considerations,  taken  together  with  what  is  known  about  the 
predictions  of  the  unaltered  direct-interaction  equations,  imply  that  the 
LHDI  equations  yield  a  Kolmogorov  inertial  range  (n  =  5/3)  and  the  associated 
Kolmogorov  dissipation  range.  Apart  from  the  appearance  of  different  time- 
correlations,  the  unaltered  direct-interaction  and  the  LHDI  equations  are 
identical;  they  have  the  same  coupling  coefficients  for  the  interaction  among 
the  wavenumber  modes.  The  unaltered  equations  give  an  inertial  range  with 


-39- 


n  •  3/2  and  local  energy-transfer  within  the  range.  It  was  shown  in  detail 
in  Ref.  1  that  this  value  for  n  is  due  to  the  appearance  of  the  Lulerian 
correlation  time  (vQk)  \  where  vQ  is  the  total  rms  velocity  component,  as 
the  effective  memory  time  for  simultaneous  triple  correlations.  It  was 
demonstrated  in  Ref.  1  that  the  Kolmogorov  spectrum  emerges  if  this  time  is 

« i 

replaced  by  (v,k) 

K 

If  the  velocity  field  has  a  Kolmogorov  inertial  range,  then  similar 
arguments  show  that  the  LriDI  scalar  equations  yield  a  5/3- law  scalar  field 
inertial  range  and  yield  Richardson's  law  for  the  relative  diffusion  of  two 
particles.  Roberts^  investigated  these  questions  for  the  unaltered  scalar 
direct-interaction  equations  and  found,  in  accord  with  the  velocity-field 
case,  that  the  Kolmogorov  and  Richardson  laws  emerged  if  the  effective 
memory  times  in  the  inertial-range  scalar  transfer  function  were  %  (v^k) 

If  the  above  arguments  are  sound,  to  whar  extent  are  the  LHDI  results 
a  justification  for  Kolmogorov's  assumptions?  Very  likelv,  the  principal 
support  they  provide  is  equally  well  sunplied  by  an  elementary  argument. 

The  LHDI  equations  can  only  demonstrate  some  self-consistenc*'  features  of 
Kolmogorov's  basic  assumptions.  But  these  features  are  already  demonstrated 
to  essentially  the  same  extent  by  the  fact  that  Kolmogorov's  dimensional 
analysis  yields  an  inertiai-range  spectrum  which  falls  within  th~  bounds  on  n 
given  in  (9.4).  A  rigorous  proof  of  the  5/3  inertial-range  law  would  involve 
much  moire  than  demonstrating  self-consistency  at  the  level  of  covariances  or 
spectra.  The  Kolmogorov  law  is  a  statistical  assertion.  If  it  is  correct, 
there  can  still  be  exceptional  flows,  or  exceptional  regions  in  a  g:  ren  flow, 
which  dc  not  obey  it.  In  order  to  estimate  the  possible  contribution  of  such 
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exceptians ,it  would  seen  necessary  to  have  some  successive  approximation 

scheme  in  which  errors  can  be  bounded. 

ifhat  the  LhDi  equations  do  provide,  assuming  again  that  they  turn  out  to 

have  well-behaved  solutions,  is  a  quantitative  embodiment  of  Kolmogorov's 

ideas,  obtained  by  precisely  defined  analytical  operations  on  the  havier- 

Stokes  equation.  It  can  be  argued  that  they  do  this  with  the  maximum 

economy  of  material  that  might  reasonably  be  expected.  The  Lagrangian 

correlation  time  seems  intimately  involved  with  Kolmogorov's  reasoning,  and 

it  is  clear  that  in  general  the  Lulerian  and  Lagrangian  correlation  times 

must  be  treated  together.  In  the  energy-range  the  two  times  are  intimately 

related.  Therefore  it  is  plausible  that  an  analytical  approximation  applicable 

to  all  spectrum  ranges  should  involve  the  full  function  U. .(x,t |r;x' ,t' |r' ). 

1  j  ~ 

Knowledge  of  the  full  function  provides  physically  interesting  information 
which  cannot  be  obtained  from  the  Euler'  1  covariance:  for  example,  the 
covariance  of  particle  acceleration. 

The  LriDI  equations  are  substantially  more  complicated  than  the  unaltered 
direct-interaction  equations  for  the  pure  Eulerian  field.  However,  they  are 
very  much  less  complicated  than  the  next-higher  Eulerian  approximation  above 
direct-interaction.  The  latter  yields  final  equations  involving  triple 
correlations  explicitly,  together  with  associated  higher  Green's  functions, 
and  consequently  involves  enormous  geometrical  complication.  Despite  this, 

it  is  only  partially  successful  in  suppressing  the  spurious  relaxation  effects 

1 

of  the  unaltered  direct-interaction  approximation. x 

If  the  details  of  the  time-correlation  functions  are  not  desired,  these 
functions  may  be  fitted  to  functional  forms  which  leave  only  certain  charac- 
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teristic  times  to  be  determined  from  the  LhDI  equations.  This  Drocedure  is 

i  *“  1 

successful  for  the  unaltered  Eulerian  direct-interaction  equations. 


10.  SGML  WEAK  POINTS  OF  THE  LriDI  APPROXIMATION 

It  was  remarked  above  that  a  model  representation  for  the  LriDI  equations 
has  not  been  found  and  that  therefore  realizable  covariance  functions  are  not 
guaranteed.  In  contrast,  the  existence  of  a  model  reoresentation  for  the 
unaltered  direct-interaction  equations  for  the  generalized  fields  pives  this 
guarantee.  The  special  solutions  so  far  found  for  tne  LHDI  equations  rive  no 
hint  of  unphvsical  behavior,  in  the  sense  of  negative  spectra  or  time- 
correlation  functions  which  do  not  fall  prooerly  to  zero  at  infinite  difference 
times,  however,  the  intricacies  of  the  equations  are  sufficient  that  further 
reassurance  would  be  comfort ine. 

The  LHDI  equations  have  been  obtained  in  this  paoer  by  an  heuristic 
modification  of  the  direct-interaction  equations  that  expr  sses  analytically 
Kolmogorov's  idea  of  examining  the  dynamics  of  the  Eulerian  field  in  coordinates 
moving  locally  with  the  flow.  Within  the  general  structural  framework  of  the 
direct -interact ion  equations,  the  alterations  of  the  triple-moment  formulas 
seem  to  be  unique.  Any  chance  in  the  rules  given  in  Sec.  6  and  8  appears  to 
result  in  less  of  one  or  more  of  the  conservation,  invariance,  and  ecui librium 
oroparties.  If  rule  (1)  of  Sec.  6  (insertion  of  S  superscripts)  is  changed, 
the  equipartition  equilibrium  solution  is  destroyed.  If  the  labeling  times 
s  are  changed  to  anything  but  t  in  anv  or  all  of  the  terms,  the  Galilean 
invariance  properties  are  lost.  However,  no  evidence  is  presented  that  the 


LhDI  approximation  is  related  to  anv  fornal  expansion  scheme  in  the  clear 
way  that  the  unaltered  direct-interaction  aoproxinat ion  is  related  to  the 
iteration  expansion  in  Dowers  of  zeroth-order  covariances  md  Green’s  functions. 
The  direct-interaction  equations  vielc  values  of  the  covariances  which,  when 
expanded,  ag ree  with  the  exact  formal  expansion  to  the  lowest  order  in  the 
zeroth-order  functions  and  which  contain  well-defined  infinite  subclasses  of 
the  higher-order  terms.  If  the  LHDI  values  are  expanded  in  the  zeroth-order 
functions,  they  also  agree  to  lowest  order,  but  so  far  no  clear  relation  of 
the  higher-order  terms  to  the  exact  formal  expansion  has  emerged.  Lven  where 
the  LiiDI  equations  give  exactly  correct  final  results,  the  structural  relation 
of  the  results  to  the  term-by-term  iteration  expansion  is  unclear.  This  may 
be  a  less  serious  lack  than  the  absence  of  a  model  representation,  because 
the  iteration  expansion,  and  all  the  other  formal  expansions  so  far  proposed 
for  the  turbulence  problem, are  almost  certainlv  divergent  and  the  collecting 
of  subclasses  of  terms  from  them  therefore  need  not  be  meaningful. 

It  should  be  noted  finally  that  the  LHDI  equations  appear  to  yield  less 
statistical  information  than  the  unaltered  direct-interaction  approximation. 

The  latter  yields  a  value  for  the  general  triple  moment  of  the  form 
^u.(x,t  |r)u.(x*  ,t’  jr’  )u  (xM,t"  ir")\ 

The  LnDI  prescription  covers  only  those  triple  moments  which  appear  in  the 
equations  of  motion  for  the  covariances;  that  is,  moments  in  which  two  of 
the  times  t,  t',  and  t"  are  equal.  If  all  three  times  are  different,  it  is 
not  clear  what  fixed  value  should  replace  the  labeling  time  s  in  carrving 
out  the  Lagrangian-history  modification. 
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APPENDIX:  EQUIPARTITION  SOLUTIONS 

Expand  the  Eulerian  field  *(x,t)  in  the  series 

*(x,t)  =  To  (t)d  (x),  (Al) 

-  *•  n  n  ~ 

where  Cr(x)  is  one  a  complete  set  of  real  orthonormal  eigenfunctions  of 
the  equation  _ 

7  ♦  +  XJn  =  °.  ( A2 ) 

n  n  n 

The  complete  set  consists  of  a  set  which  vanishes  everywhere  on  the  boundaries 

and  a  set  whose  normal  derivative  vanishes  everywhere  on  the  boundaries.  Let 

n  {=  1,2,...)  label  the  eigenfunctions  in  order  of  nondecreasing  1^. 

Eq.  (2.9)  with  «c  =  0  gives 

$  (t)  =  T  A  (t)*(t),  (A3) 

n  ltbl  nn.  a 

where  the  A  (t)  are  functionals  of  u(x,t|t).  The  sum 
nm  —  — 

[^(t)]2  =  /0(x,t)J2d3x  (t.  u) 

is  a  constant  of  motion  for  any  initial  condition.  Suppose  that  only  and 

»  are  excited  initially.  Constancy  of  (A4)  at  the  initial  instant  then  shows 
in 

that  the  A^^t)  identically  satisfy 

A  (t)  =  -  A  ( t ) ,  (AS) 

nm  an 

whence  the  system  satisfies  the  Liouville  theorem 

(t)/3*  (t)  =  0.  ( A6) 

It  follows  that  the  equipartition  probability  density 

£>(♦)  *  exp[-(  const  J^*2]  (A7) 

is  an  absolute  equilibrium  distribution.  (A7)  is  to  be  understood  in  the 

following  sense.  All  *  and  t>  for  n  or  re  greater  than  N  (N  as  large  as  desired) 

n  m 

A 

are  dropped  from  (A3)  and  (A4),  yielding  a  finite  system.  If  G^^tjt')  is 
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14 

the  response  matrix  of  the  i^’s,  then  it  follows  that 

(*  (tlMt'l)  =(G  (t;t')>  (t  >  t’),  (A8) 

\  n  in  /  x  nm  ' 

where  the  average  is  over  the  suitably  normalized  distribution  (A7),  Trans¬ 
formation  back  to  x  space,  followed  by  the  limit  N  -*•  •  then  gives 

¥(x;t it;x' ,t' |t’ )  =  G(x;t|t;x’ ,t* |t* )  (t  >  t'),  (A9) 

and  (5.15)  then  follows  from  the  fact  that  there  if  no  r  or  r'  dependence 
when  <  *  0. 

The  inviscid  equipartition  solution  for  the  Eulerian  velocity  field 
satisfying  (3.21)  in  an  arbitrary  domain  may  be  established  by  a  similar 
procedure,  using  an  appropriate  set  of  eigenfunctions.  [Note  that  this 
solution  contains  both  shear  and  compressive  velocity  compcnents  but  that 
the  compressive  compont  its  are  static  and  do  not  affect  the  dynamics  of  the 
shear  components.]  The  technique  of  Ref.  14  new  establishes  (8.15)  for 
r  •?  t,  r'  *  t*.  The  relation  may  then  be  established  for  general  r  and  r* 
by  applying  the  technique  of  Ref.  14  to  (3.18).  This  readily  shows  that 
for  any  r,  the  field  u(x,tjr),  considered  as  a  function  of  t,  is  in  eaui- 
partition  equilibrium  and  satisfies  a  fluctuation-relaxation  relation  like 
that  of  the  Eulerian  field. 
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FIGURE  CAPTIONS 

Fig.  1.  Evolution  of  u(x,t|r)  in  the  (t,r)  plane  from  the  initial  time  tQ. 
The  point  (t,r)  may  be  either  above  or  below  the  diagonal. 

Fig.  2.  Propagation  route  in  the  (t,r)  plane  of  the  perturbation  described 
by  G(x,t |r;x' ,t' jr' ).  The  points  (t,r)  and  (t',r')  may  be  indeDendentlv 
either  above  or  below  the  diagonal,  but  r  >  r*  always. 

Fig.  3.  Paths  of  evolution  and  propagation  associated  with  the  direct- 
interaction  contributions  to  S  and  H.  The  points  (t,r)  and  (t’,r*)  may  be 
independently  either  above  or  below  the  diagonal.  In  Fig.  3a,  r  >  r’  or 
r  <  r*.  In  Fig.  3b,  r  >  r'. 

Fig.  4.  (a)  Paths  of  evolution  associated  with  the  direct-interaction  or 

LHDI  contribution  to  S.  and  M.  .  (b)  Paths  of  propagation  associated  with 

in  in  *  r 

the  contribution  to  H.  and  N.  . 

m  in 
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FIG.  I 
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(b) 

FIG.  4 


